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We monitor the relaxation ol internal stresses in a fractal colloidal gel on cessation of flow and 
find a weak power law decay, (Ji ~ t~ a over 5 decades of time where a w 0.07. The system exhibits 
physical aging of the elastic modulus, G' ~ t 13 , with f3 « a. Imposition of zero stress after waiting 
time t w results in strain recovery as the system relaxes without constraint. Remarkably, recoveries 
at different t w can be shifted to construct a master curve where data are scaled vertically by l/ai(t w ) 
and plotted horizontally as (t — t w )/tt^ where fj, w 1.25, indicative of a super-aging response. 

PACS numbers: 82.70.Dd,82.70.Gg,83.80.Hj,61.43.Hv 
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Many out of equilibrium systems display a slow evolu- 
tion of their dynamics and properties in time, towards an 
eventual stationary or equilibrium state. This process is 
referred to as "aging" , and it is the hallmark of many dis- 
ordered materials, ranging from molecular jpolymer and 
spin glasses, to colloidal gels and glasses [l|, |2j, & U, [5( . 
Although the origin of the aging response varies widely 
in these systems, they share the commonality of a slow 
recovery of some canonical structural or dynamic quan- 
tity after a departure from equilibrium, usually initiated 
by a rapid change in a corresponding control parame- 
ter. For example, after a sudden temperature quench in 
polymer glasses, evolution of free volume and heat capac- 
ity lend themselves to both conceptual and experimental 
investigation and have proved highly successful as met- 
rics for the recovery behavior of polymers Q. In spin 
glasses, below the spin glass transition temperature T c , 
removal of the external magnetic field leads to a slow 
and non-exponential decay of remanent magnetization 
on long time scales Q- In colloidal systems, the pic- 
ture is somewhat murkier Q . Attempts to correlate the 
evolution in dynamics to measurable changes in struc- 
ture have so far been unsuccessful The concept of 
internal strains, and the resultant stresses, due to out of 
equilibrium locations of constituent matter is a potential 
starting point, for systems with soft potentials. It has 
been advanced with some success as a model for the dy- 
namics of metallic glasses [l(| wherein the distribution of 
stresses in the system defines the energy landscape that is 
traversed in the slow return to equilibrium during aging. 

In this Letter, we consider strain recovery as a possible 
structural metric for the aging of an attractive colloidal 
system and examine the relationship between strain re- 
covery and internal stress present in the system. We 
study a dispersion which exhibits rapid gelation via the 
formation of a fractal network on cessation of shear flow 
and follow the evolution of the shear modulus in time 
after gelation. We separately conduct experiments to 
monitor the relaxation of internal stresses after flow ces- 



sation, as well as the strain recovery of the system after 
various waiting times during the stress relaxation. We 
find that the strain recovery behavior as a function of 
system age can be collapsed into a single curve using a 
vertical shift factor that is inversely proportional to the 
internal stress at the start of the strain recovery and a 
horizontal shift that scales inversely with elapsed time. 

Our system consists of dilute dispersions, 2-6 wt.%, 
of carbon black particles (Cabot Vulcan XC72R) in a 3:1 
mixture of mineral oil and tetradecane ( Aldrich Chemical 
Co.). Samples are prepared by dispersing the particles in 
the solvent under vigourous vortexer mixing for 2 min- 
utes, followed by sonication for at least 30 minutes. Sam- 
ples are studied in a stress controlled rheometer (Anton- 
Paar MCR301) using a 50 mm 1° cone-plate geometry. 
The instrument is mounted on an air table (Newport) 
to eliminate ambient mechanical vibrations which could 
perturb the sample. Two separate sets of measurements 
were performed. The first was of the time evolution of 
the system modulus. The second characterized internal 
stress relaxation and subsequent strain recovery. The 
protocol for the modulus evolution consists of 4 steps: 
(1) a high shear rate rejuvenation step at 7=1000 s _1 
which completely erases the flow history of the material 
and ensures a reproducible starting point [llj, (2) pre- 
shear at the shear rate of interest, 7=100 s _1 , for 1200 
s to achieve a steady-state viscosity, (3) cessation of flow 
(the shear rate was reduced from 100 s _1 to in 0.1 s) 
and (4) measurement of the modulus as a function of 
time at u> = 1 rad/s and 7 = 0.1%, within the linear 
viscoelastic regime. We set t — at the end of step 3, 
where flow stoppage leads to rapid system gelation. The 
protocol for the second set of experiments consists of five 
steps. Steps 1-3 are identical to that of the first protocol. 
In step 4, the sample is maintained under a quiescent or 
zero strain rate condition (the system is stationary). The 
finite, decaying stress required to maintain the station- 
ary condition is the internal stress [Hj]. We measured 
the internal stress after cessation of flow as a function of 
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time for a waiting period t w . In step 5, at t = t w , the 
zero strain rate is replaced by a zero stress condition, al- 
lowing the system to undergo strain recovery, analogous 
to that in traditional creep recovery experiments. The 
recovered strain is measured as a function of time via 
the displacement of the rheometer tool. Measurements 
consisting of steps 1-5 were conducted for a range of t w 
from 1 to 10 4 s. The fourth step of our protocol provides 
a waiting time during which the sample ages quiescently. 
In typical rheological and dynamic scattering investiga- 
tions, system properties are not monitored during the 
waiting step. The protocol applied here, however, per- 
mits correlation of the system properties measured after 
the waiting time with the internal stress dynamics char- 
acterized therein. 

Under the conditions of the experiment, the carbon 
black particles interact via an attractive Van der Waals 
potential that results in the formation of a col- 
loidal gel at rest. Gelation in these systems is extremely 
fast [12|, ll4| such that a substantial elastic modulus, well 
in excess of the viscous modulus, is observed at the earli- 
est measurable times after cessation of flow, ~ 0.1 s. The 
elastic modulus displays a slow power-law increase with 
time, Figure [H where G' ~ with (3 « 0.06. This slow 
growth is due to thermally driven structural evolution of 
the fractal gel, as observed for analogous systems [1 51 ] . 
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preserved across the numerous iterations of the stress re- 
laxation and strain recovery measurements, as shown for 
a 4 wt.% sample in Figure [3] Similar consistent behavior 
was observed for the other concentrations underlining the 
efficacy of the rejuvenating step in effectively eradicating 
the flow history of the sample. 
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FIG. 2: Relaxation of internal stress as a function of time t 
across different samples studied. 



-i — i 1 1 mii| — i — i 1 1 1 ■ ■ ■ | — i — i 1 1 1 ■ 1 1 1 — i — i 1 1 1 ■ 1 1 1 — i — i 1 1 1 ■ 1 1 1 




ul_ 



ul_ 



ul_ 



0.1 



10 100 
t(s) 



1000 



10000 



FIG. 1: Evolution of elastic modulus as measured via oscilla- 
tory shear with 7=0.1% and lu=1 rad/s. 

Conversely, the internal stresses that are established in 
the system during the rapid gelation on cessation of flow 
exhibit a weak decay in time. The data are well fit by 
a power law, where the stress decays as <7j ~ t~ a over 
5 decades of time with a « 0.07, as shown in Figured] 
The close correspondence between the scaling exponents 
of the stress relaxation and the aging of the gel modulus 
suggests that the relaxation of internal stresses may be a 
key indicator for the aging of the system. Strain recov- 
eries were recorded after various stress relaxation dura- 
tions, t w . The trajectories of the stress decays are well 



FIG. 3: Internal stress decay during different waiting times 
t = t w for the 4 wt.% sample. The stress is monitored up 
to the different values of t w shown in the legend in separate 
iterations of the experiment. 

The imposition of zero stress after waiting times t w of 
1, 10, 100, 1E3 and 1E4 s resulted in strain recoveries 
ranging from 1 — 5%, measured over 2000s, Figured] At 
very short times, oscillations are present in the data due 
to inertially driven ringing of the sample. Such "creep 
ringing" is commonly observed in elastic gels when sub- 
jected to a stress impulse [lf3|. Here, the impulse origi- 
nates from the internal stress of the sample itself which 
drives an initial fast recoil of the tool that is quickly 
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damped out by the viscosity of the sample. Thereafter, 
the material undergoes a slow strain recovery as struc- 
tural rearrangements occur during aging. Younger sam- 
ples, those at smaller t w , display an initial rapid strain 
recovery which appears to asymptote towards a long time 
plateau value. By contrast older samples exhibit slower 
initial recovery from a short time plateau, followed by a 
gradual rise at long time. Samples of intermediate age 
show an inflection in their recovery. Remarkably, strain 
recovery data from different sample ages can be shifted 
to construct a single master curve at each composition, 
Figure [5] 
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FIG. 4: Strain recovery after pre-shear at 100 s~ for different 
samples after different waiting times t w . The concentration is 
mentioned in the figure. The times mentioned in the legend 
are different values of t w . 

The vertical shift factors are simply inversely propor- 
tional to the value of the internal stress at the start of 
the strain recovery, b(t w ) — l/(Xi(t w ). The horizontal 
shift factor a(t w ) — where the best overlay is ob- 
tained for [i = 5/4. The vertical shift of the data by the 
inverse of the internal stress is strikingly simple. It is rep- 
resentative of an elastic strain that is recovered quickly 
at short times upon the start of the measurement. At 
long times, there is a far slower and continuous recovery 
of strain that persists out to the longest times measured, 



100 e — i 1 1 1 1 1 1 1 1 1 1 — a 




10"" 10' 7 10" 6 10" 5 10"* 10" 3 10" 2 10"' 10° 10 1 10 2 10 3 10 4 
(t-t w )a(t w ) (s) 

FIG. 5: Strain recovery master curves assembled from data 
taken as a function of t w . Squares, circles, triangles, diamonds 
and stars denote 2 wt.%, 3 wt.%, 4 wt.%, 5 wt.% and 6 wt.% 
concentrations respectively. a(t w ) = t^, 5 ^ 4 , b(t w ) — (Ji(f„) _1 . 
The times mentioned in the legend are different values of t w . 



over 1000 s. We have ruled out external influences such 
as tool drift or room vibration in the system. It should be 
noted that the complex modulus and thus the dynamic 
viscosities of the samples are large, > 200 Pa for the 
data of Figure 21 so any small perturbations from exter- 
nal sources would be quickly dissipated. The shifting of 
dynamical responses for older samples to short rescaled 
time is commonly encountered in soft glassy systems 
0, H3] and has been related phenomenologically 
to the elapsed time rescaling first successfully advanced 
for polymer melts plj . In these systems, however, either 
simple aging (n = 1) or sub-aging (n < 1) are observed, 
whereas here we obtain the best overlap for /i > 1, in- 
dicative of a super-aging response. The master curves in- 
dicate that the dynamics of strain recovery at long times 
should be understood from the behavior of young materi- 
als, and vice-versa, that the dynamics at short times are 
exhibited within the observation window by older sam- 
ples. Such a display is counter intuitive if viewed simply 
in the context of recovered strain as marking the prox- 
imity to long-time equilibrium. In that sense, older sam- 
ples would show more asymptotic tendency than younger 
samples and not the other way around as observed. This 
display, however, depends on the duration of the waiting 
time relative to the timescale for complete relaxation. 
For short t w , we can view the response purely as a re- 
flection of the slow vs. fast dynamics of old vs. young 
samples. Thus we expect elapsed time superposition to 
hold as it indeed does. The near uniform scaling of the 
internal stress relaxation with time across the different 
compositions studied (Figure El (a) = 0.07 ± 0.002) en- 
ables the strain recovery to be simply scaled in terms of 
waiting time alone where the vertical shift factor is now 
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b <~ The data of Figure [5l across 5 compositions, can 
all be rescaled with the same shift factors with very good 
fidelity. 

Internal stresses are known to persist for long times in 
a variety of disordered systems 2^, 23, 24, 25|. Stochas- 
tic, local stress relaxation events have been invoked as a 
concept to explain the unusual intermittcncy and ballistic 
dynamics observed in various soft materials 22, 24, 25j |. 
The origin of these stresses is in the non-equilibrium lo- 
cation of system constituents. In colloidal systems with 
soft potentials, the sudden nature of an ergodicity break- 
ing transition results in the arrest of particles away from 
preferred locations, that is those where the gradient of 
the inter-particle potential is minimized. As a result, 
forces exist among particles locally, and are propagated 
throughout the system along particle contact chains or 
branches of the fractal network in the case of gels. The 
departure of the system from equilibrium is encoded by 
the distribution of local displacements of particles with 
respect to their equilibrium positions. Aging occurs via 
structural rearrangements that minimize these internal 
stresses as particles conduct a thermally driven explo- 
ration of their local potential energy landscape. In this 
framework, it is clear that the evolution of local strain, 
and in response, local stress, chart the aging process. In 
a gel which is formed by the slow aggregation of freely 
diffusing species, such local displacements and stresses 
are directionally random and the resultant macroscopi- 
cally observable stress is zero, or vanishingly small. In 
the present system, the application of shear implies that 
the displacement of particles from their preferred equilib- 
ria and the deformation of particle clusters will be biased 
along the flow direction. The rapid gelation that occurs 
on cessation of flow then results in residual stresses that 
are macroscopically non-zero. As expected, the stress 
acts counter to the direction of the deformation, so the 
sign of the stress is reversed if the direction of the shear 
flow is changed 12j. 

The correspondence between the relaxation of internal 
stress and the increase of the gel modulus point strongly 
to the identification of the internal stress state as an in- 
dicator for the aging response of the system. From a 
microscopic perspective, due to the deforming effect of 
the shear flow before gelation, structural rearrangements 
involved in aging are anisotropic and thus give rise to 
macroscopically observable strain recovery. The dynam- 
ics of this strain recovery then serve as a well rationalized 
metric of the aging behavior of the system, as demon- 
strated here. At long times, the total recovered strain, 
[y(t— t w )]t =0 oi should be smaller for older samples, which 
have undergone longer stress relaxation under zero strain. 
The master curves however imply no significant difference 
in the long-time asymptotic values of recovered strain for 
different sample ages. This is an implication of the small- 
ness of the experimental time scales relative to that for 
the full aging of the system. This point is underscored 



by the continuous power-law dependence of the internal 
stress relaxation over 5 decades of time. 



The authors gratefully acknowledge L. Cipelletti, V. 
Trappe and D. Weitz for helpful discussions and NSF 
funding via CBET-0828905. 



[1 

[2 
[3. 

K 

[5, 

ir 

[§ 
[9 
[10 

[11 
[12. 

[13 

[14; 
[is; 

[16 
[17 

[is; 

[19 

[20 

[21 
[22 

[23 

[24 

[25 



* ajay.negi@yale.edu 
1 chinedum.osuji@yale.edu 

] W. Kob and J. L. Barrat, Phys. Rev. Lett. 78, 4581 

(1997) . 

] C. Angell, K. Ngai, G. McKenna, P. McMillan, and 

S. Martin, J. Appl. Phys. 88, 3113 (2000). 
] S. Hilgenfeldt, S. A. Koehler, and H. A. Stone, Phys. 

Rev. Lett. 86, 4704 (2001). 
] L. Cipelletti and L. Ramos, Curr. Opin. Colloid Interface 

Sci. 7, 228 (2002). 
■] L. Cipelletti, L. Ramos, S. Manley, E. Pitard, D. Weitz, 
E. Pashkovski, and M. Johansson, Faraday Discuss. 123, 
237 (2003). 

] J. D. Ferry, Viscoelastic Properties of Polymers (Wiley, 

New York, 1980). 
] K. Binder and A. P. Young, Rev. Mod. Phys. 58, 801 
(1986). 

■] G. B. McKenna, T. Narita, and F. Lequeux, J. Rheol. 

53, 489 (2009). 
'] G. C. Cianci, R. E. Courtland, and E. R. Weeks, Solid 

State Commun. 139, 599 (2006). 
! ] T. Egami, S. J. Poon, Z. Zhang, and V. Keppens, Phys. 

Rev. B 76, 024203 (2007). 
] A. S. Negi and C. O. Osuji, Rheol. Acta (2009). 
] C. O. Osuji, C. Kim, and D. A. Weitz, Phys. Rev. E 77, 

060402 (R) (2008). 
] V. Trappe and D. A. Weitz, Phys. Rev. Lett. 85, 449 

(2000) . 

] S. R. Raghavan, H. J. Walls, and S. A. Khan, Langmuir 

16, 7920 (2000). 
] L. Cipelletti, S. Manley, R. C. Ball, and D. A. Weitz, 

Phys. Rev. Lett. 84, 2275 (2000). 
] C. Baravian and D. Quemada, Rheol. Acta 37, 223 

(1998) . 

'] M. Cloitre, R. Borrega, and L. Leibler, Phys. Rev. Lett. 

85, 4819 (2000). 
,] X. F. Shi, A. Mandanici, and G. B. McKenna, J. Chem. 

Phys. 123, 174507 (2005). 
'] G. Ovarlez and P. Coussot, Phys. Rev. E 76, 011406 
(2007). 

i] G. R. K. Reddy and Y. M. Joshi, J. Appl. Phys. 104, 

094901 (2008). 
] L. C. E. Struik, Polymer 38, 4053 (1997). 
:] L. Ramos and L. Cipelletti, Phys. Rev. Lett. 87, 245503 

(2001) . 

.] R. R. Hartley and R. P. Behringer, Nature 421, 928 
(2003). 

J M. Bellour, A. Knaebel, J. L. Harden, F. Lequeux, and 

J.-P. Munch, Phys. Rev. E 67, 031405 (2003). 
] R. Bandyopadhyay, D. Liang, H. Yardimci, D. A. Ses- 
soms, M. A. Borthwick, S. G. J. Mochrie, J. L. Harden, 
and R. L. Leheny, Phys. Rev. Lett. 93, 228302 (2004). 



